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Abstract 

The non-abelian Einstein-Born-Infeld-Dilaton theory, which rules the dynamics of 
tensor-scalar gravitation coupled to a su(2)-valued gauge field ruled by Born-Infeld 
lagrangian, is studied in a cosmological framework. The microscopic energy exchange 
between the gauge field and the dilaton which results from a non-universality of the 
coupling to gravity modifies the usual behaviour of tensor-scalar theories coupled to 
matter fluids. General cosmological evolutions are derived for different couplings to 
gravitation and a comparison to universal coupling is highlighted. Evidences of cos- 
mic acceleration are presented when the evolution is interpreted in the Jordan physical 
frame of a matter respecting the weak equivalence principle. The importance for the 
mechanism of cosmic acceleration of the dynamics of the Born-Infeld gauge field, the 
attraction role of the matter fluid and the non-universality of the gravitational cou- 
plings is briefly outlined. 



1 Introduction 



In the description of the very beginning of the universe, well before the Big 
Bang nucleosynthesis, field theoretical models are to be considered instead of 
the usual hydrodynamical description of matter. Those kinds of models, inspired 
by high-energy physics, have lead to numerous progress in modern cosmology, 
trying to solve various problems from cosmic acceleration or the flatness problem 
to the magnitude of the cosmological constant or the existence of topological de- 
fects. The motivation of this paper lies in two questions, amongst many others, 
that raise from the description of the very first moments by high energy physics. 

First is the question whether large scale massless gauge fields can play any 
interesting role in cosmology. Indeed, such fields could have existed in the early 
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universe before the phase transitions of spontaneous symmetry breaking but, if 
they were ruled by usual Yang-Mills (YM) conformally invariant dynamics, their 
primeval excitations have probably been swept away by inflation. This point 
motivated some authors to study the impact of the Born-Infeld (BI) type mod- 
ification of gauge dynamics, suggested by string theory, on cosmology (see pQ). 
The BI lagrangian breaks down the scale invariance of the gauge fields beyond 
some critical energy, and therefore it is not obvious to conclude directly on the 
becoming of such gauge fields during and after an inflation period. Furthermore, 
it was proved in P that such gauge fields of Bl-type cannot provide any cosmic 
acceleration on their own although they can mimic a fluid of negative pressure. 
Before going further on this question, it is therefore of first importance to study 
more deeply in a cosmological context the interaction between gauge and scalar 
fields as suggested in models inspired by high-energy physics. 
The second question is to see what happens to a possible scalar sector of grav- 
itation during the cosmological evolution. Indeed, string theories predict the 
existence of the dilaton a Lorentz scalar partner to the tensor Einstein 
graviton as low-energy limit of bosonic actions. This large theoretical frame- 
work provides a physical background for tensor-scalar modification [3j of gen- 
eral relativity, in which gravitation is mediated by long-range scalar field acting 
in complement of the usual spin 2 gravity fields. Although this question has 
been widely studied when dilaton - or more generally tensor-scalar theories - 
is coupled to matter during radiation and matter dominated era, the case of a 
microscopic field model which would not be coupled universally to gravity, as 
suggested in string theory, has been less considered. In particular, how does the 
interaction between scalar and gauge fields modify their respective dynamics 
and the resulting cosmological evolution will be the main subject of the present 
paper. 

But before going any further, let us locate the present work in the existing lit- 
terature. In this paper, we will focus on cosmological solutions of the Einstein- 
Born-Infcld-Dilaton (EBID) equations for flat spacetimes. Cosmologies with 
large scale massless homogeneous and isotropic gauge fields with gauge group 
SU{2) and ruled by usual YM dynamics have been studied for a long time 
( 01 HI |Sj ) • The gravitational instability of flat spacetimes filled with such 
gauge fields was studied in Generalisation to higher gauge groups have also 
been studied ^01 \U\ in the case of flat and closed cosmologies. The Einstein- 
Born-Infeld cosmology with non-abclian gauge fields deriving from gauge group 
SU{2) has been studied thoroughly in PQ| for flat, closed and open spacetimes 
for any value of the cosmological constant. The minimal coupling of large scale 
cosmological gauge fields and scalar multiplets has been studied in ^1 
The Einstein- Yang-Mills-Dilaton (EYMD) equations for flat cosmologies and a 
special case of non-universal coupling to gravity, have been derived in ^] ^j] 
where the authors highlighted the energy exchange between the dilaton and the 
gauge fields and briefly discussed its effects on inflation, entropy crisis and the 
Polonyi problem (domination of a nearly massless dilaton at late times). How- 
ever, they did not propose a complete solution to the EYMD field equations 
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in the cosmological framework which could allow to address completely these 
issues. Moreover, they did not discuss the general influence of non-universal cou- 
pling to gravity as well. The EYMD system was also studied in J7| in the case 
of closed Friedmann-Lemaitre-Robertson- Walker (FLRW) with a static gauge 
field 1 and vanishing dilaton potential and cosmological constant. 

As we told before, tensor-scalar cosmologies have been widely studied, with 
a large spectrum of applications for physical cosmology: inflation, primordial 
nucleosynthesis, Cosmic Microwave Background, ... The question of the con- 
vergence of tensor-scalar theories to general relativity during the cosmological 
evolution has been widely studied in ^3 ^3 1201 El and references therein. For 
the so-called "Einstein conformed frame" , where the gravitational and scalar 
fields have pure spin 2 and spin dynamics respectively, the scalar sector of 
gravitation disappears naturally during cosmic expansion due to its coupling to 
matter. The cosmological evolution of the dilaton emerging from string theory 
has been studied in |22j . 

In this paper, we will consider the cosmological evolution of the dilaton coupled 
directly to a large-scale non-abelian gauge field ruled by BI dynamics which go 
beyond the scale invariance of YM theory. A non-universal coupling to gravity, 
as suggested in preceding works, will lead to quite different results to the usual 
coupling of tensor-scalar theories to a fluid. For example, when the gauge fields 
are governed by YM scale-invariant dynamics, the scalar sector of gravitation 
remains directly coupled to the gauge fields although they mimic a radiation 
fluid. However, it is well-known that tensor-scalar theories decouple from radi- 
ation (except during phase transition). Through both numerical computations 
and analytical solutions, we will show how the dilaton evolution is modified by 
non- universal coupling to the metric. This will lead to remarkable consequences 
for cosmology. 

The structure of this paper will be as follows. In section 2, we establish the 
general field equations for the EBID cosmology. In section 3, we first remind 
the reader about BI cosmology as was studied in [J. The non-abelian BI cos- 
mology can be split in two extreme regimes depending on the energy density of 
the gauge field compared with a critical scale introduced in the BI theory. For 
large energies, the gauge field is shown to mimic a fluid of negative pressure with 
p/p = —1/3 while in the low-energy limit the scale-invariant YM dynamics is 
retrieved and the gauge field looks like a radiation fluid. On the other hand, we 
also remind the reader about dilaton cosmology, studied in ^J|22]. As the equa- 
tions are to be established and solved with pure spin degrees of freedom in the 
Einstein frame, the interpretation in terms of the Jordan physical frame is also 
recalled. The case of universal coupling is solved there with usual properties of 
tensor-scalar theories. In section 4, we will focus on the strong field limit of the 

1 Due to this particular topology, the EYMD system does not reduce to pure Einstein- 
Dilaton field equations when the gauge field is static. 
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EBID system with general coupling to gravity and in section 5 the low-energy 
limit which consists of a generalised version of the EYMD system appearing in 
|14l ITS] is treated. In section 6, we analyse a complete cosmological evolution 
of the EBID system and present evidences for possible cosmic acceleration in 
the physical Jordan frame. This frame is defined with respect to a pressureless 
matter fluid that has been added to the gauge sector. The acceleration is shown 
to resist to the attraction provided by matter and appears to be intrinsically 
related to the non-universality of the coupling to gravity. Finally, we conclude 
in section 7 by some perspectives to the present work. 



2 Field Equations of Einstein-Born-lnfeld-Dilaton Cosmology 

Most of the interest of field models, for example deriving from string theory 
in the low energy- limit, comes from their non-universal coupling to the gravity 
fields g^ v and its scalar counterpart <j) (each type of matter field has in general its 
own coupling function to the dilaton, see [221 )■ This results in a violation of the 
weak equivalence principle and therefore the gravitational interaction of these 
microscopic field models is different from an usual tensor-scalar theory where 
the weak equivalence principle is usually assumed. Without imposing such a 
violation of the weak equivalence principle at a microscopic scale, field models 
would not be different than considering a tensor-scalar theory in presence of 
a fluid with the equation of state of the considered fields. Therefore, we will 
make use of a general form of the action for the non-abelian Einstcin-Born- 
Infeld-Dilaton system, that takes into account a possible violation of the weak 
equivalence principle. This action writes down 

S = J {-^B -^cpd^-V(0)-A 4 (<f > )£ BI (B 2 (<f>)g^,A^^gd i x 

+S m [ip m ,C 2 (<t>) 9fiU ]- (1) 

In this action, the gravitational interaction is described by the scalar curvature R 
and the dilaton cf>, k being the "bare" gravitational coupling constant and A(<p), 
B{<f>) and C(4>) being three different coupling functions of the dilaton to matter. 
The first two illustrate the coupling of the gauge sector to the volume form and 
to the Einstein metric g^ v inside the Born-Infeld lagrangian C bi (where A^ are 
the non-abelian gauge potentials) and the last coupling function C(<p) is related 
to another type of matter ruled by the action S m . Another parametrisations 
of Einstein-Born-lnfeld-Dilaton action were considered in the litterature, for 
example with A = 1 and B = cxp(fc/2</») in [H El HI El and with B = 
A 2 = cxp(fc/20) in j2HJ • The non-abelian gauge interaction is described by the 
Born-Infeld lagrangian build upon the field strength tensor 
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(with T a are the generators of the gauge group under consideration 2 ) and its 
dual tensor F^ v . Indeed, this lagrangian, denoted by Cbi in ||TJ, is defined as 

Cbi = e c (H-l) 

where e c is the Born-Infeld critical energy and £?(</>) is the dilaton coupling 
function which is equal to e k l 2 ^ when non-perturbative effects are not taken 
into account (see ^]E3E2])- I n this case, k will be called the dilaton coupling 
constant. Throughout this paper, we will assume the Planck system of units, in 
which h = c = 1 and G = nip 2 , with the Planck mass mpi = 1.2211 x 10 19 GeV 
and the gravitational coupling constant is k — 8irG. We have also set the gauge 
coupling constant to unity, as it actually defines a system of units for the dilaton 
field </>, provided the dilaton being massless (V — 0). The Born-Infeld critical 
energy e c defines the scale above which non local effects of string theory arise 
and where the scale invariance of the gauge fields is broken. For example, in the 
low-energy limit e c — > oo of the Born-Infeld part of the action we recover 
the usual, conformally invariant, Yang-Mills (YM) lagrangian density for the 
non-abelian gauge field: 

V~gC YM = -V^A 4 (0)B- 4 (0)iF^Fr • (3) 

In this paper, we will focus on cosmology and therefore we will adopt the pre- 
scriptions of the cosmological principle which states that the spatial sections of 
our Universe are homogeneous and isotropic. For the sake of simplicity, we will 
also restrict ourselves to the case of flat space-times, which constitutes how- 
ever a very nice approximation of the present universe and its early stages as 
well. The metric describing such spacetimes is the one of Friedmann-Lemaitre- 
Robertson- Walker (FLRW): 

ds 2 = ~N 2 {t)dt 2 + a 2 {t) (dr 2 + r 2 d0 2 + r 2 sin 2 6dip 2 ) (4) 

where a(t) is the scale factor and N is the so-called lapse function of the hamil- 
tonian Arnowitt-Deser-Misner approach to general relativity. This function can 
be fixed by a specific choice of time coordinate (gravitational gauge freedom). 
The symmetries prescribed by the cosmological principle we assume impose that 
the dilaton scalar field <p depends only on time. 

A remarkable fact about non-abelian gauge fields is that they admit non-trivial 
homogeneous and isotropic configurations at the opposite of their abelian U{1) 
counterparts (see ^ |S] and references therein for a complete discussion of the 
SU (2) case). The main reason for that is because only the gauge invariant quan- 
tities such as the field strength tensor have to exhibit the symmetries explicitly, 

2 Gauge indices will be noted as bold latin letters. 
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while the gauge potentials can be symmetric up to a gauge transformation (see 
|2"H 0] for more general gauge groups). As a result, the energy can be dis- 
tributed amongst the different gauge degrees of freedom while the stress-energy 
tensor remain compatible with the maximal symmetry of the space-time back- 
ground. In this paper, we will restrict ourselves to the case of su(2)— valued 
gauge potentials, for which the ansatz 

A = A^dx" = a{t)T m dx m (5) 

of the connexion one-form A makes the gauge invariant quantities satisfying the 
required symmetry (see [5] ) . The remaining dynamical degrees of freedom of 
the gauge potential are now expressed by the field cr(t). However, our results 
will not depend on this particular choice as the ansatz above can be generalised 
to higher gauge groups (see ^HE])- In the equation above, the generators T m 
of the Lie algebra of the gauge group SU (2) are to be taken in the coordinate 
dependent basis of the gauge degrees of freedom space as follows: 



T r = sin 9 cos ipTi + sin 9 sin ipT-z + cos 9T$ 
09 J 



Te = ±T, 



T - - T 



sin 9dtp 



with Tj = jdi the usual basis of the Lie algebra su(2) (<7j being the Pauli ma- 
trices) with the following standard normalization conditions and commutation 
relations: 

tr {T a T h ) = - J ab , [Ta,T b ] = ie c ah T c 

The ansatz (JSJ is of course independent of the particular choice of the action 
for the gauge field, as was shown in p^. 



The symmetries implied by the cosmological principle therefore allow us to 
write down Q as an effective one dimensional action, after integrating over R 3 
and dividing by the infinite volume of its orbits: 

Seff = Jdt^^ + ^-V^Na'-Nah^^ill-l^+S^K (6) 

where a dot denotes a derivative with respect to the time t and where 1Z is given 

« - lA-s^f-^-^-o^*"^. CD 

V e c \a?N 2 a 4 J e 2 c a 6 N 2 v ; 

where a if) is the gauge potential as defined in (JSJ. Following 1 , it is also 
convenient to write TZ as 

K = v 7 ! - tVTTA (8) 
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with 



r = 



3d 2 £- 4 (0) 



A = 



3a 4 B" 4 (0) 



(9) 



e c a 2 N 2 

From the action JfjJ) and relation J7J , it is straightforward to write down the field 
equations for the Einstein-Born-lnfeld-Dilaton system by varying this action 
over the following degrees of freedom: N, a, <p and a. First, the Euler-Lagrange 
equation for the variable N gives the hamiltonian constraint 



^ + V{<j>) + e c A\^) {V-l) + Pn 



(10) 



which we will refer to as the Friedmann equation. In the previous equation, p* 
stands for the energy density of the matter fluid ruled by S m and the function 
V is defined in terms of T and A in @ as 



V = 



i-r 



(ii) 



The careful reader should have noticed that, after varying over N, we set the 
gravitational gauge to N = 1, meaning that we work with the synchronous time 
coordinate (another convenient choice for the study of the gauge dynamics in 
the Yang-Mills regime is the conformal gauge N = a as it naturally exhibits the 
conformal invariance). The Friedmann equation allows us to define the Born- 
Infeld effective energy density of the gauge field as a generalisation of what was 
proposed in P 

PBX = (V-l). (12) 

The Euler-Lagrange equation for the scale factor a gives the acceleration equa- 
tion: 



(v(<t>) 



1 



4?) + e c A\cp) (V- 1 - 1) - - ( Pm + 3 Pm ) 



(13) 



where p m stands for the pressure of the additional matter fluid. This allow us 
to define the Born-Infcld effective pressure 



Pbi 



-A 4 (<P) (3-V-2V- 1 ) 



(14) 



and the equation of state 



Pbi 
Pbi 



Pbi 



Pbi 



(15) 



for the gauge part of the EBID system as in pQ. Here we see that the su(2)~ valued 
gauge fields ruled by Born-Infcld lagrangian can be represented by a fluid with 
an equation of state that varies continuously from — | , when the BI energy den- 
sity is much larger that the 'critical field' A~ a (4>)pbi 3> e c , to i at low energies 
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A~ 4 ((j))pBi <C e c . These two extreme regimes correspond to a gas of Nambu- 
Goto strings in three spatial dimensions on one hand (strong field limit) and 
radiations on the other (weak field limit). The transition between these regimes 
occurs at vanishing pressure when the BI energy density is of order of the BI 
critical energy scale e c . At low energies (e c — > oo), the gauge field behaves like 
radiation as expected because the BI lagrangian reduces to the conformally in- 
variant Yang-Mills one. 

It is also important to notice that there is no cosmic acceleration with the 
metric as long as the dilaton is massless. Indeed, the highest value of a that 
can be achieved in this frame is identically zero fsee !13|) . in the limit of the pure 
Einstein-Born-Infeld system at high energies (0 = 0, Xbi = — However, we 
will see that cosmic acceleration may appear once we examine the behaviour in 
another frame. 

By varying the action @ with respect to the dilaton </>, we find the Klein- 
Gordon equation: 

a acp 

-l(4>) (Pm ~ 3p m ) (16) 

where a(0) = ^J^, 0(0) = and 7 (0) = The key point of 

the physics in the EBID system lies in the fact that the dilaton field couples 
differently to the gauge sector of the theory (last term) depending on the values 
of the coupling functions A and B. Although the cosmological dynamics of the 
gauge fields ruled by Born-Infeld lagrangian can be regarded as a fluid with an 
equation of state given by l|15|l . the coupled dynamics of the dilaton and the 
gauge field does not reduce in general to a scalar-tensor theory with this fluid 
as background. This is only the case when we have a universal coupling to the 
metric g^ v , i.e. when A = B. In the general case, there exists a non-trivial 
energy exchange between the dilaton and the gauge sectors of the theory that 
will dominate at late epochs as we shall see further. 

Finally, the Euler-Lagrange equation for the gauge field a gives the Born-Infeld 
equation that rules the gauge potential dynamics: 



<T n a ( 1 > 

\—V- 2 ~2-cj \--V- 2 )+2(f>& 
a z a \ 2 



2a(0)|-/3(0) (2^ + l-V- 2 



= 0- 



(17) 

This equation is essentially the same that in ^ except from the coupling term 
proportional to <p which accounts for the direct energy exchange at the micro- 
scopic level between the fields. It is important to notice that we did not assume 
any direct coupling between the gauge field and the additional matter fluid, 
which will allow to treat them separately. 
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Indeed, from equation Q17)). and following pp, it is possible to derive an en- 
ergy conservation equation for the BI density : 

Phi = -2- PBI PBI t^l? + MtifecA'M (ft 1) 
a Pbi + e c ^l 4 (0) 

-2p(^e c A\4>) (^±1 -2V + 2ft) ■ (18) 

Now that we have derived the complete set of the EBID field equations for 
cosmology, we propose the reader to briefly review some basic features of Born- 
lnfeld cosmology on one hand and dilaton cosmology on the other. In the 
rest of this paper, we will only consider a massless dilaton (i.e., vanishing self- 
interaction potential V((f>) = 0). 

3 Born-lnfeld and Dilaton Cosmologies 
3.1 Non-Abelian Born-lnfeld Cosmology 

The Non-Abelian Born-lnfeld cosmology in various spacetimes with different 
values of the curvature and the cosmological constant was described in detail in 
1 . The field equations governing these models are those of the previous section 
with a vanishing dilaton <p = = 0, constant coupling functions A = B = 1 and 
no additional matter fluid p m = 0. The equation (|18|) for BI energy conservation 
can be written as 

a p B i + 2e c 

PBI = -2- Pbi ; (19) 

a pbi + e c 

which admits a first integral: 

a 4 pbi (pbi + 2e c ) = C (20) 

where C is a positive constant. In the strong field limit, psi 3> e c , the BI energy 
density rcdshifts as psi ~ a ■ while in the weak limit, pbi <C e c , we retrieve 
the radiation behaviour pbi ~ & characteristic of the conformal invariance of 
the gauge field at such energies. This allows to treat separately the spacetime 
evolution and the dynamics of the gauge field. Although the complete analyti- 
cal solutions for both gravitational and gauge sector were derived in [Q, let us 
illustrate simply the main features of this cosmological model. 

First, the strong field limit psi ^> e c corresponds to V ^> 1. In this limit, 
the acceleration equation (Tf^|l reduces to 



a k 
a 6 



whose general solution is 



a(t) = a* sin ( \j 
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where we set a(0) = and a* is the value of the scale factor at the time 
t* = \ — \ (in Planck units). Therefore, the cosmic expansion starts with a 



zero acceleration at the singularity. Then, setting V 3> 1 in equation l|17fl brings 

a 

(7 — — (7 = 

a 

which shows that a scales as a. Near the singularity, the behaviour of the gauge 
field is therefore 

and the gauge potential a starts at rest. 

Then, in the weak field regime, psi <C e c and V ~ 1 (e c — > oo). This limit 
corresponds to the Einstein- Yang-Mills cosmological solution studied in The 
conformal invariance of the gauge field in that regime yields that the scale factor 
behaves like in the radiation-dominated era : a(t) ss ^/t (in synhronous time). 
On the other hand, the energy conservation equation i|18fl now reduces to 



2 • 2 , 

a a + a = 



4 __C_ 
3e c 



which can be integrated in terms of the Jacobi elliptic function. Moving to the 
conformal time coordinate dt = adrj, we find 



r(r)) = £ 1/4 cn (£ x/ \ 



-1 



where cn(u, k) is the Jacobi elliptic function and £ = In synchronous time, 
the gauge potential a oscillates with a fixed amplitude and a growing period. 

More generally, it is possible to derive a general solution for the gauge potential. 
Let us rewrite the first integral i12()ll in terms of V as 



p =f + i?- (21) 

Using the definitions and © , the previous equation may be integrated 
to give the gauge potential (in the conformal gauge dt — adr/): 

= al\fe~ c cn {V^-q\ -l) • 

Figure^illustrates the evolution of the scale factor, the gauge potential and the 
equation of state during the expansion of a non-abelian Born-lnfeld universe. 
The figures correspond to the numerical integration of equations and l|17fl 
with <p = <j) = and A = B = 1. During numerical evolution, we monitor the 
violation of the hamiltonian constraint l|l(J|) (see the appendix for more details 
on integrating the EBID system). In this case, this violation does not exceed a 
part on 10 -12 . Initial conditions at U = were assumed such as pBi{ti) = 100e c , 
a(ti) = a, = 1 and &(ti) = (e c = 10~ 4 x m^). 
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Fig. 1: Illustration of Non-Abelian Born-lnfeld cosmology : a) scale factor b) 
gauge potential c) equation of state of the Born-lnfeld "fluid" 



3.2 Dilaton Cosmology 

Dilaton cosmology can retrieved from our fundamental equations by setting 
a equal to 0. Indeed, tensor-scalar theories can be written in the so-called 
"Einstein" conformal frame: 

S= hJ dix ^9~*i R * - 2<?r^M + S m [tp m , C 2 {^)g%] , (22) 

where k is therefore the "bare" gravitational coupling constant, S m is the action 
for the matter fields ip m , tp — y/ n/2<f) and <?*„ the "Einstein" metric tensor 
which corresponds to basic gravitational variables with pure spin 2 propagation 
modes. This metric is measured by using purely gravitational rods and clocks 
and allows to account for the dynamics in a simpler way 3 than a observable 
frame in which the metric tensor g^ v is universally coupled to matter fields ip m . 
This frame is called the " Jordan- Fierz" frame in which the action 122(1 can be 
written 

S =\J d 4 x^UR- C ^r"d^d v ^+S m [i; m ,g liV }, (23) 

where R is the curvature scalar build upon the physical metric g^y which is 
measured using non-gravitational rods and clocks and where u> (<E>) is called 
the coupling function. The scalar field $ now gives the effective gravitational 
coupling constant. In this frame, matter fields evolves in the same way that they 
could do in general relativity because the action of matter does not depend 

3 In particular, this frame represents gravitation with its pure scalar and tensor degrees of 
freedom and the limit of general relativity is not singular. 
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explicitly on the scalar field <f>, as matter couples universally to the physical 
metric g^. Einstein and Jordan frames can be linked together through the 
conformal transformation 

g ta , = C 2 {<p)g* pa , (24) 

and the following relation between the scalar field <$ in the Einstein frame (pure 
spin dynamics) and its counterpart $ in the Jordan physical frame: 

$ 1 = kC 2 fo>) , (25) 

with ip — n/2(j>. The energy density and pressure of the matter fluid in both 
frames are related by 

pm = C 4 p™ (26) 
Pm = C 4 p~ m , (27) 

where p m and p m represent these quantities expressed in the Einstein frame. At 
the opposite of scalar-tensor theories, we did not assume here an universal cou- 
pling to the metric tensor g^. Indeed, the action (Tjp) reduces to a tensor-scalar 
theory in the Einstein frame |23) only when A{4>) = B{4>) = C(<j>), i.e. when 
the weak equivalence principle applies. Let us now remind the reader about 
the behaviour of such tensor-scalar theories in the presence of a background 
cosmological fluid. When the matter fields ip m are represented under the ap- 
proximation of a perfect fluid, the dynamics of the scalar field is ruled by (see 

[El EH EH E2): 

2 <fi" + (1 - A) iff + (1 - 3A) 7 fo>) - 0, (28) 



(3 

where 7(y>) = d '"J^ and A = p/p is the equation of state for the cosmological 
fluid. In the previous equation, a prime denotes the derivative with respect to 
the variable p = \n(a/ai). The action of the cosmological fluid is thus to damp 
the dynamics of the scalar field while it is rolling down some effective potential 
depending on the coupling function. Furthermore, the scalar field now has an 
effective, velocity-dependent, mass of 

mfo) = (29) 

where the field has a limiting speed cp 1 < v3 for which its effective mass diverges. 
This relativistic limit corresponds to the case where the energy density of the 
background fluid is negligible compared to the kinetic energy of the scalar field 
(the universe is dominated by the kinetic energy of the scalar field). 

3.3 The universal coupling for EBID cosmology 

Let us now turn back to the EBID system we wrote in the previous section 
and focus on the gauge sector only by setting S m = 0. If we now assume a 
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universal coupling to the metric by setting A = B, we now have for the 
dilaton equation l|16|) : 

4> + 3^0 + a{4>) {pbi - ipBi) = 0, (30) 

where 

Pbi - 3p BI = 2e c A 4 {cp) (V + V- 1 - 2) . 

In this case of universal coupling, we recognize the equation for the scalar field in 
presence of a background fluid for general tensor-scalar theories . The equations 
1101) , (|13l) and 11GI) can be simply solved in the Jordan frame by using the results 
on the non-abelian Born-lnfeld cosmology (previous paragraph). The equation 
<|17|) for the gauge field in the Einstein frame can also be written 

u + 2^V- 2 -2^&(l- V- 2 ) - 2a{<t>)j>a (l - V' 2 ) = • (31) 
a z a \2 J 

Therefore, in the weak energy regime psi <C e c (V ~ 1), the gauge field under- 
goes a conformally invariant dynamics (non-abelian radiation, psi — 3pbi) and 
decouples from the scalar field. 

In a radiation-dominated universe, where the equation of state is A = p/ p = 
1/3, the dynamics of the scalar field is given by the following solution to l|28|l 

(cf. [HH3I2I]): 



ip(p) =cp oc -v^ln Ke- p + Vl + K 2 e- 2 P , (32) 

where the integration constant K is determined from the initial velocity ip'(p — 
0)=ip' a : 



K 



'2 



This should correspond to the low-energy limit of the Born-lnfeld field equa- 
tions when a universal coupling is assumed: the scalar field velocity in p-time 
should be damped to zero by the cosmological expansion. It should be noticed 
that when there is no universal coupling A ^ B, we do keep an energy exchange 
between the dilaton and the gauge fields and the usual dynamics of tensor-scalar 
theories will be modified. Moreover, when there is no universal coupling, the 
energy exchange between the gauge potentials and the dilaton field will prevent 
the dynamics to be purely dictated by the solution ()32fl. As we shall see in sec- 
tion 5, this solution will accurately describe the early epochs of evolution when 
the field is almost relativistic. However, at late times, energy transfer between 
dilaton and gauge fields will substantially alter the dynamics. 

Once again, let us assume a universal coupling and consider the strong field 

4 We also have pbi = A^psi (see !12l for the relation between the energy density expressed 
in the Einstein and Jordan frames (quantities with a ~). 
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limit of the Born-lnfeld system where we have A = —1/3. The dilaton equation 
<|28|) now becomes 



9 



+ - V ' + a(ip) = • 



3 -ip' 2 3' 

Let us now write down for the dilaton coupling function 



A\4>) = 

= 

a((p) = 
|3 + 2w($)| = 



2k - 3/c 2 
2k 2 



'2k 



(33) 

(34) 
(35) 

(36) 
(37) 



and face the simplest tensor-scalar theory of Brans-Dicke type (tp = J k/24>). 
Using CONVODE [35], it is possible to find an analytic solution for tp' under 
the following implicit form when we use the coupling function (|34|l : 



(6^ 2 -8)(p + p ) = \/3^1n 



-V3 



V3 



- 2 In 



'2 



(38) 



where A = k/\/2K and po some integration constant. When p 
attractor for ip', namely 



(3A + 2(p'YJ 

oo, there is an 



<P'(P - 



oo) = — 
' 2 



3 k 



'2k 



(39) 



Therefore, there is also a maximum value for the dilatonic coupling constant 
k for which the attractor corresponds to the relativistic limit for the dilaton 

kmax = ' (40) 

In the non-relativistic limit ^ \/3, equation l|33|l now becomes 

p" + 2<p' + 3A = 0, 
which can be solved easily to give 



3 k 
2 



2 k 



P 



(41) 



where we assumed tp i = 0. We can see that, due to the constant potential term 
in equation (|33|l . the value of the dilaton field goes to — oo (+oo if k < 0) with 
the time- variable p. However, the gauge energy density will also decrease with 
time and finally the assumption of strong field will be no longer true as pbi 
becomes less than the critical energy e c . At the end of the evolution, we should 
retrieve the radiation case for which the solution in case of universal coupling 
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was described above. Once again, a non-universal coupling to the metric yields 
modification of these behaviours as we shall sec further. 



Now that we have recalled the main features of Born-Infeld and dilaton cos- 
mologies as well as EBID system with universal coupling, let us now discuss 
how the non universal coupling in the general EBID system will modify a tensor- 
scalar cosmological picture. This will be done in three steps: in the following 
section, we will focus on the strong field limit where the gauge field mimics a 
Nambu-Goto string gas ; then on the low-energy limit which corresponds to the 
Yang-Mills regime for the gauge fields and finally to the general cosmological 
evolution where transition between both regimes occurs. 



4 The Strong Field Regime 

As we have seen earlier, the strong field limit is reached when the BI critical 
energy e c can be neglected with regards to the gauge field energy density. Setting 
V > 1 into the EBID field equations (JTUJl, 02J) and (JTHJl and (O with p m = 
p m = and V = 0, we find: 



+ 



cf> 2 + e c A\$) 



3-0 + 2e c A 4 (0)7> [-/3(<j>) 



2a((f>)] = 



-O- + 20O- [2a(<f>) - 3/3(0)] = 0- 



(42) 

(43) 
(44) 
(45) 



where we used ^ = 1 — T and assumed T <C 1 which will be verified afterwards 
by the agreement between analytical and numerical solutions. 



In the following, we will assume the exponential coupling function Ij34(l for the 
sake of simplicity. However, the qualitative analysis will be valid for any cou- 
pling functions. 



In terms of the new time variable 



p = In I — 



(where defines the initial zero value for p) we can combine equations B42JI to 
1(31} to obtain 

P "~ + \<p' + (-/%) + M<P)) - 0, (46) 



3-tp' 2 3 1 



4 The Strong Field Regime 



16 



where <p = y/ ' K/2tfi- In order to particularise, we can set now a = (A = 1) and 
B(ip) = exp(fc/v / 2Ktp) and find for the dilaton equation 

«" ■ 2 - ' =0 (47) 



3-^' 2 3^ ^ 

which admits a solution similar to the case of a tensor-scalar theory with Nambu- 
Goto string gas (l33l) : 



(6.A 2 - 8) (p + po) = -VSAln 



21nf y,a " 3 a V (48) 



where .4 = kj\rlK- When p — > oo, the attractor for is now exactly the 
opposite of the universal coupling case: 

¥/(p-oo) = §-j|=. (49) 

The maximum value for the dilatonic coupling constant k is the same as before 
(equation (|4(J|) ~) and the non-relativistic limit can be obtained from (|41|l with 
an opposite sign. The constant potential term in equation l|46|) is negative so 
that the value of the dilaton field is pushed toward +00 (—00 if k < 0) with 
the time-variable p, as long as the gauge field remains in the strong field limit 
Pbi > e c . 

Figure El a) illustrates the evolution of the dilaton velocity with respect to 
the time variable p in the strong field limit. The trajectory has been com- 
puted numerically (solid line in Figure [21 a)) from the integration of the full 
EBID system with the initial conditions indicated in the caption (with A = 1 
and B(ip) = exp(k/\/2K(p)). Also shown is the analytical approximation of the 
strong field limit given by equation (|4*8t represented by big dots. 

Figure [21 b) gives the evolution of the dilaton field along the cosmic expansion 
for the case A — 1 and B((p) — exp(k/ y/2nip) . Starting with a negative velocity, 
the dilaton is damped to a minimum before being accelerated to infinite values 
(with k > 0). Fortunately, as the BI energy density will decrease with time, 
the strong field limit psi ^> £c will soon be no longer valid. We shall see fur- 
ther that, in the Yang-Mills limit, the same coupling to gauge fields will bring 
the dilaton to infinitely negative values. Of course, similar conclusion can be 
found when k is negative. Figure |3 c) gives the behaviour of the scale factor 
in the case discussed here. As the dilaton field is relativistic at infinitely- low 
times p — > —00 and therefore dominates the energy content of the universe, the 
expansion starts with an infinite rate, breaking the Ll renormalisatiorC that was 
done in simple BI cosmologies. Finally, let us focus on the gauge sector of the 
EBID system. In the strong field limit, the features of Born-Infeld cosmologies 
for the gauge field are conserved: it starts at rest, damped by the cosmic ex- 
pansion, before entering the oscillation regime of the weak energy Yang-Mills 
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Fig. 2: Evolution of the EBID system in the strong-field regime with the p-time 
variable a) <p f b) dilaton field (j) c) scale factor a(t) d) Energy density of 
the BI gauge field (ai — 1, pBi{ai)/e c = 10 10 , e c = 10~ 4 x rrip^ k = 5, 
$ = -1.73, (5i7/i7 < 10" 13 ) 
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limit. Figure [21 d) illustrates the evolution of V for the numerical solution (solid 
line) with dilaton compared to the evolution in a simple BI universe with same 
initial BI energy density (dash-dotted line). This holds for the particular cou- 
pling A — 1 and B(ip) — exp(fc/\/2Ktp). If we now transpose A and B, it is 
easy to see from l[46l) that the attracting value for ip' will be twice the value of 
the universal coupling (|39fl . Therefore, when the coupling to the volume form 
is weaker than the coupling to the metric 5 , the scalar field behaves just the 
opposite way than in the universal coupling (with A = —1/3). Up to this point, 
we have obtained both analytical and numerical solutions for the strong field 
limit, pbi S> e c , of the EBID system. We have also explained qualitatively the 
effects of non-universal coupling to the Einstein metric on the dynamics of 
the scalar field. Let us now turn to the weak field regime in which the gauge 
sector is ruled by Yang-Mills lagrangian. 



5 The Weak Field Regime: solutions of the 
Einstein-Yang-Mills-Dilaton system 

The weak field regime of the EBID system is reached when the BI energy density 
A A (ip)pBi becomes much smaller than the critical energy e c . In this case, the 
lagrangian ruling the gauge sector takes the usual Yang-Mills form J3J), which 
gives for the spatially homogoneous and isotropic gauge potentials Jsj: 



C 



YM 



N 2 af 



(50) 



Taking into account the limit e c — > oo (and V ss 1) into the EBID field equations 
EH, lO, OH and D3J (with V = 0), we obtain: 



a 
a 



K 

3 

K 

~3 

& 

a 



-A\<f>)B-\ct>) 



°-A\4>)B-\& 



+ 2- 



QA^)B-\<p) (afo)-/3(0) 
-& + 44& (a((f>) - f3{4>)) = 0- 



(51) 

(52) 
(53) 
(54) 



These equations constitute the Einstein-Yang-Mills-Dilaton system, and the spe- 
cial case of A = 1 and B(<f>) = exp(fc/20) can be found in [141 115) . In this paper, 
the authors highlighted the importance of the energy exchange between the dila- 
ton and the Yang-Mills field. Indeed, this coupling yields a new force term in the 
field equation for the dilaton 153[) and the gauge field l[54l) which disappear in 
case of universal coupling {A — B and a — (3). With non-universal coupling the 



For example, in the case A = 1 and B(ip) = exp(k/v2Kip) we just discussed 
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gravitation is now sensitive to the force term when it is coupled to Yang- Mills 
radiation although its equation of state is those of radiation (see (J2HJ)- Let us 
now describe the dynamics of the different fields in the EBID system for this 
low-energy regime. 

First, let us move to the p-time variable p — ln(a/ai) and use the accelera- 
tion and Friedmann equations l|52[). i|51[l to rewrite (|53|l as 



+ ^--2{a{4>)-m) -22,4 =°> ( 55 ) 



Z-tp' 2 3 v ' & 2 a 2 , .. 

with 

The energy exchange term in & 2 a 2 — a inside relations (|53[l and (|55|l is in gen- 
eral oscillating, due to the self-coupling of the non-abelian gauge field (term in 
2a 3 /a 2 in H54|l). A way to handle this easily is to replace it by an effective source 
term which would account for the average effect of the gauge field oscillations. 
Let us therefore proceed to the following replacement 

a a 



f « N (56) 

with H some constant expressing the effectiveness of the energy exchange be- 
tween dilaton and gauge fields in the weak field regime. This constant H can 
for instance be estimated numerically by computing the average of the driving 
term in (|55|l over one period. Equation (|55|l is the same field equation than l|28|l 
for the tensor-scalar theory of the dilaton but now with a non-vanishing force 
term due to our averaging of the gauge oscillations. By averaging the gauge 
oscillations, we obtain a similar equation to the strong energy limit (equation 
l|46|) ') seen in the previous section. Therefore, we can use the same procedure as 
before : if we set A = 1 and B{4>) = exp(fc/2</>), we can propose the following 
implicit solution for <p'\ 



~2){p + p a ) = -V? l A\n 



V3 



9' 



V3 



-Inp^V), (57) 



where A = 2Nfc/\/2K. Once again, the p-time derivative of the dilaton field ip' 
evolves towards the following attractor 

If (p — > OOJ = — D 



. 2k 

and the maximum value allowed for the dilatonic coupling constant k for which 
the dilaton remains relativistic (ip'^ — > — v3) is 
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It is important to notice the opposite sign between the attractors of the strong 
and weak field regimes which will have important consequences on the cosmo- 
logical evolution of the dilaton. In the non-relativistic limit, <p 2 < 3, we find 
the following solution for the dilaton 

<p= -3A(e- p +p-l) +<fi (58) 

and we see that the dilaton tends to — oo, if the dilatonic coupling constant is 
positive (N > 0). 

Figure compares our analytical solution i|57f) coming from an averaging ap- 
proximation (dashed line) to a numerical solution of the full EBID system (solid 
line). In the YM regime, the velocity of the dilaton cp' appears oscillating, 
around average values given by the approximation (|57|) with H = 1/3. This 
value numerically appeared to account for the average behaviour of the EBID 
system in the low-energy limit for a wide range of parameters (k, e c , tpi or ip'A. 
Therefore, averaging the oscillations of the source term to about a third of their 
amplitude seems in very agreement with numerical solutions. Therefore, the 
attractor for the p-time derivative of the dilaton is now 

^(p - oo) = -2-* (59) 
while the maximum value allowed for the dilatonic coupling constant is 




When A — 1 and B((f) = exp(fc/2</>), the p-time derivative of the dilaton ap- 
pears to converge to a constant negative value which is directly given by the 
non-relativistic approximation l|58|) which is valid when k is small compared 
with y/ln) (dotted line). When B = 1 and A(c/>) — exp(fc/20), it is obvious 
from equation l|55|) that the attractor has exactly the opposite value. Therefore, 
the case of non-universal coupling A ^ B is quite different to what happens 
in a usual tensor-scalar theory: in a radiation-dominated universe, the p-time 
velocity of the dilaton freezes to zero (see the solution l|32|l). Indeed, this will 
make the dilaton field diverging after an infinite amout of time as can be seen 
in Figure In this figure, we represented the evolution of the dilaton in the 
numerical solution with A = 1 and B(<fr) = exp(fc/2</>) (solid line) to the solu- 
tion (|32|l for the same field in a radiation-dominated universe. More precisely, 
the dilaton tends to — oo. As a conclusion, although the gauge field looks like 
radiation at a large-scale level, the non-universal microscopic coupling between 
dilaton and gauge sectors finally dominates at large redshift and freezes the en- 
ergy contribution of both sectors in such a way that none of these components 
completely dominates. Figure[S]presents the evolution of the gauge field veloc- 
ity & with the scale factor. The gauge field appear to be damped by its coupling 
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Fig. 4: Evolution of the dilaton field as a function of the p-time variable in 
the weak field limit (Same parameters as above) 
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Fig. 5: Evolution of the gauge field velocity as a function of the scale factor in 
the weak field limit (Same parameters as above) 

to the dilaton and in fact the whole gauge sector looses energy at a rate fixed 
by (JUS). 

Figure El shows the departure of the scale factor from the radiation solution: 



Q>rad 

for the numerical solution presented in this section. We see that the departure 
is important when the scalar field is dominating at early times (t < 5 x 10 tpi) 
then finally converges to a slightly less strong expansion at late times (a « 
0.99 x cirad), when the equilibrium has been reached. 

Before going further, let us summarise the cosmological evolution of the 
fields constituting the EBID system in the YM regime. First, the dilaton is 
damped until its velocity is attracted to a negative (resp. positive) value for the 
particular coupling A — 1 and B(<f) = exp(fc/2</>) and k > (resp. B = 1 and 
A(4>) = exp(fe/2<jf>)). However, it should have been damped to rest if it would 
have been plunged in a radiation-dominated universe with universal coupling. 
Because its velocity has been attracted to a negative value, the dilaton field 
will eventually diverge linearly to — oo (k > 0). This is exactly the opposite 
situation of the strong field limit that was presented before where the coupling 
term drives the dilaton to infinitely high values. In a general situation where the 
gauge field starts with an energy much higher than the BI critical energy and 
then cools down to YM dynamics, one should expect that the dilaton reaches 
some extremum value (tp' — 0) during the transition. This will be treated in 
the next section. 
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Fig. 6: Evolution of the departure from a radiation-dominated universe of the 
weak field limit during expansion (Same parameters as above) 

6 General Cosmological Evolution 

Let us follow in detail some typical cosmological evolutions of the EBID system 
for various couplings. Starting at singularity, the dilaton is in general relativis- 
tic (|y| — > oo and ip' 2 — ► 3). The expansion therefore begins at a = with an 
infinite rate and the gauge field dynamics is dominated by the non-local effects 
induced by the BI non-linearity (pbi £c)- As energy is exchanged between 
the dilaton and the gauge field in the strong-field regime, the dilaton velocity 
is attracted to some value depending on the coupling functions, as illustrated 
in figure During this phase, the velocity of the dilaton in p-time is indeed 
a positive (negative) constant when A = 1 and B((f>) = exp(fc/20) (B = 1 and 
A(4>) = exp(k/2(j>) or the universal coupling A = B). When the gauge energy 
density has decreased to the BI critical energy (pbi ~ e c ), the dilaton veloc- 
ity leaves the strong-field attractor to enter the YM low-energy regime. The 
epoch of this transition varies according to the coupling functions considered 
(see Figure [7J . It then moves to the low-energy attractor by accomplishing 
damped oscillations around the analytical solutions proposed in the previous 
section (with A ^ B) or is damped to vanishing velocities when there is univer- 
sal coupling (see section 3). With non-universal coupling (A ^ B), the value 
of the dilaton field reaches some extremum (<// = 0) on its way to the second 
attractor. This extremum is unavoidable as we have seen that the attractors of 
the dilaton velocity in the strong and weak field regimes are of opposite signs 
and its velocity will therefore vanish at some time during the transition between 
these two attracting regimes. It is also interesting to examine the evolution of 
the gauge field energy density. Figure |H1 represents the gauge field energy den- 
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Fig. 7: Evolution of the dilaton velocity with p-time for three different couplings: 
(3 = and a = k/2 (solid line) ; a = and (3 = k/2 (dashed-dot) and 
universal coupling {A = B, dotted line) (oj = 1, pBi(cii)/e c = 10 4 , 
e c = 1(T 4 x m A Pl , fc = 2, $ = 0.1, < lO" 10 ) 




Fig. 8: Evolution of the physical energy density p with p = ln(exp(fc/20)a) for 
the curves in Figure [7| (Same parameters as above) 
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sity (V — 1) related to the curves in Figure [7| With universal coupling A = B, 
we retrieve a cosmological evolution given by equation l|19|) and its first integral 
<|2U|) : pbi ~ aT 2 in the strong field regime and Pbi ~ a -4 at low energies. The 
assumption of non-universal coupling now leads to different evolutions of the 
gauge field energy density which in fact are given by the more general energy 
conservation equation l|18|) • The differences between the evolutions comes from 
the particular trajectories illustrated in Figured 

Now that we have reviewed the main features of a general full cosmological 
evolution from high to low energy regimes, it is also important to discuss the 
evolution of the observed cosmological parameters like the scale factor, the Hub- 
ble expansion rate or the accelerating parameter in the Jordan physical frame. 
In order to define such a frame, we will introduce an additional matter pressure- 
less fluid which will verify the weak equivalence principle. The energy density 
p m of this fluid in the Einstein frame is related to the physical energy density 
p m through the relation (|26[) . The coupling function C((f>) to ordinary matter 
defines now our 11 observable" Jordan frame by 

= C 2 (<j>)g^ ■ (61) 

In the Jordan frame obtained by the previous conformal transformation, the 
energy density p m of the matter fluid is ruled by the same conservation laws as 
in general relativity. This is true because there is no direct interaction between 
the gauge sector and the additional matter fluid and therefore they decouple 
from each other. The field equation for the gauge potential a does not need to 
be modified as we do not assume any direct coupling with the pressureless fluid. 
We will now consider the field equations for the EBID system we have written 
in section 2 for the Einstein metric g^ v with a pressureless fluid p m = 0. 

The observable scale factor in the Jordan frame will be given by 

a = C(<t>)a, (62) 
while the synchronous time in the Jordan frame is denoted by 

dt = C(4>)dt ■ (63) 
Then, the Hubble expansion rate can be derived directly 

H = ^- i = C- 1 (4>)(H + y(cf>)<j>), (64) 

where H = a/ a is the Hubble parameter in the Einstein frame, a prime de- 
noting a derivative with respect to p = ln(o/aj) (Einstein frame) and j((j>) = 
dhxC(<f>)/dcf>. The acceleration parameter q in the Jordan frame can be written 



a a 
a 



a ^^^0 2 a + 7(0)0a + 7(^)0a + ^ (0 + 07(^)0) ^ (65) 
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where a dot over a quantity expressed in the Jordan frame means a derivative 
with respect to the synchronous time t in that frame. Although there is no 
possibility of a cosmic acceleration (q > 0) in the Einstein frame (unless one 
considers a non- vanishing potential, see relation 1(13(1 1. this does not rule out 
a possible acceleration for the observable scale factor given by J?HJ|. Indeed, 
the existence of fluid (constituted by our Born-Infeld non-abelian gauge field) 
that violates the weak equivalence principle will result in a possibility of cosmic 
acceleration. To show that this is actually the case even in presence of matter 
fluid which would make the tensor-scalar theory converging to general relativity 
if taken alone, we will use the coupling function 7(1/5) = 92. Figure El illustrates 
evolutions of the acceleration parameter in the Jordan frame q(p) , given by 1(65(1 
as a function of the Jordan scale factor a. The solutions presented here are 
expanding universes (H > 0). Four different couplings have been considered, 
including the simple case of A = B ^ C which corresponds to different couplings 
of gravitation to the gauge and the matter sector. The pressureless fluid energy 
density at start has been chosen to dominate the BI energy density by more than 
one order of magnitude. As the gauge field starts in the strong field regime, its 
energy density will scale approximately with a -2 , depending on the coupling 
functions (see also Figure |SJ|. Therefore, the gauge sector rapidly dominates 
the energy content of the universe. The dynamics of the dilaton is as described 
earlier: after having moved to the strong field regime attractor, the transition to 
YM dynamics occurs and the dilaton quickly moves to the low-energy attractor. 
During this transition, acceleration appears in the Jordan frame defined by the 
pressureless fluid as indicated in Figure El We see also that any violation of the 
weak equivalence principle (here by taking C ^ A, C ^ B or C ^ A = B) lead 
to cosmic acceleration even with tensor-scalar theory that would alone converge 
to general relativity. An important condition for cosmic acceleration is to have 
a repulsive force term in the dilaton equation. Therefore, the EBID system with 
a non-universal coupling to gravitation offers the interesting possibility to build 
a scenario for dark energy or inflation. For the curves represented here, the 
acceleration periods are shorter than in a usual LCDM or quintessence model 6 . 
and therefore a more complete study should be done to determine if it is possible 
to explain distance-redshifts measurements with EBID fields. It should also be 
noticed that an EBID dark energy scenario would predict a finite period of 
acceleration. Indeed, as the gauge field will recover a YM dynamics at the 
end of its evolution, its energy density will finally scales as a~ A and will finally 

6 To give an idea on how these universes accelerate, we remind the reader about the following 
values of the acceleration parameter for various energy content: 

qlradiation) = —1 
q(relativistic <j>) = —2 
9(A) = 1 
q(ghost) = 2 

where A tends for the cosmological constant (the de Sitter solution to which inflation is usually 
matched as an exponential expansion). 
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Fig. 9: Evolution of acceleration parameter in the Jordan frame for four different 
couplings: j3 = and a = k/2 (solid line) ; a = and (3 — k/2 (dashed- 
dot) ; universal coupling (a = (3 = 7 = ip, dotted line) and a = (3 = kef), 
7 = ip (dashed-dashed, k < 0) (dj = 1, PBi{o-i)/^c — 10 5 , e c = 10~ 4 x 
m% „ fc - 4, Jff/ff < 10- 10 , ^ = 0.1, p m (ai)/pB/(ai) w 100) 



be dominated by a pressureless fluid. Therefore, the questions whether cosmic 
acceleration (q > 0) can occur, with which intensity and for how long seems 
to depend on both initial energy distribution, the value of the dilaton coupling 
constant and the critical BI energy scale e c . More work should focus on that 
point to sec if it would be possible to build a physically relevant quintessence 
model with the EBID field equations. However, the perspectives of cosmic 
acceleration in the Jordan frame do exist and this looks particularly interesting 
for our view of modern cosmology. 

7 Conclusion 

The non-abelian Einstein-Born- infcld-Dilaton model provides an interesting frame- 
work, motivated by string theory, to study the impact of large scale non abelian 
gauge fields on tensor-scalar theories of the gravitational interaction. In this pa- 
per, we focused on the cosmological evolution of an homogeneous and isotropic 
configuration of these fields in a flat background. The microscopic coupling 
between the dilaton and gauge fields induced by non-universal coupling to the 
metric leads to an energy exchange between both gravitational and gauge sec- 
tors that will alter the usual dynamics of tensor-scalar theories (for which the 
coupling is universal). 
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As in non-abelian Born-Infeld cosmology, we considered two different regimes 
depending on the gauge energy density compared to the critical energy that 
parametrizes the Born-Infeld lagrangian. We have derived both analytical and 
numerical solutions to describe the cosmological evolution of the whole system. 
In the case of non-universal coupling, the gravitational scalar field no more de- 
pends on the equation of state of the gauge field and the dynamics is altered as 
follows. 

In the particular case of a Brans-Dicke theory, in which non-perturbative terms 
for the dilaton are not considered, we have shown that the energy exchange 
resulting from the particular couplings damps the dilaton to a frozen non- 
vanishing velocity. In the high energy regime of the gauge dynamics, the at- 
tracting value for the velocity is positive when the gauge field couples more 
to the metric than to the volume form. The opposite situation happens in 
the low-energy regime where the gauge field is ruled by Yang-Mills lagrangian. 
Therefore, in a general cosmological evolution where the gauge field cools down 
to low-energies, there is a transition between the two attractors. Their values 
are directly proportional to the value of the dilatonic coupling constant. 

However, it is well-known from experimental tests of the gravitational theory, 
especially the determination of the post-newtonian parameter 7, that the value 
of the coupling ui n is at least of order 500, roughly 10~ 3 for an (see f° r 
a recent estimation of the post-newtonian parameter 7). The influence of the 
dilaton potential is also important to consider. Furthermore, the constraints on 
the weak equivalence principle obtained by the tests on the universality of free 
fall exclude a violation of this principle that would exceed a part in 10~ 12 . One 
can therefore argue that the effect of such non-universal couplings should be 
neglected. But, if the violation of the weak equivalence principle only applies 
to large-scale fields which do not couple to ordinary matter and whose distri- 
bution is roughly homogeneous, their energy density on our scales is far beyond 
experimental reach and the violation could be hard to exhibit. 

The interesting possibility introduced by such a violation is a cosmic accel- 
eration in the physical frame associated to ordinary matter. In this work, we 
build a first simple model based on our treatment of the EBID field equations 
that exhibits periods of accelerations in presence of ordinary matter verifying 
the weak equivalence principle. The acceleration has been shown to resist to 
the attracting property of the accompanying matter and seems to be a general 
feature of a non-universal coupling to gravitation. Furthermore, this model re- 
spects the weak energy condition p + 3p > in the frame of the physical degrees 
of freedom (the Einstein frame). The Born-Infeld dynamics of the gauge field 
plays a crucial role in this kind of dark energy model by ensuring a late aris- 
ing of this mechanism (when the gauge field mimics a Nambu-Goto string gas) 
and even predicts an end to the dark energy domination (when the gauge field 
looks like radiation). An interesting perspective to this work would be to use 
this remarkable feature of non-abelian Born-Infeld gauge fields to build physical 
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models for quintessence and, maybe, inflation. 

In conclusion, we can say that our study of non-abelian Born-Infeld gauge fields 
coupled to to tensor-scalar gravity opens new and interesting perspectives for 
the question of the attraction to general relativity as well as other crucial topics 
of modern cosmology such as inflation or dark energy. 
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Appendix : numerical integration of the EBID system 

Here we give some details of the numerical integration of the full EBID system 
we use to illustrate this paper. In order to integrate the system of equations 
, (fTUfl . ijTT>|) and (fr7|) . we choose the following procedure. First, we rewrite 
equations (|13[) . (|16[1 and l|17|l as a system of six first-order ODE's and keep the 
hamiltonian constraint (|10|l to check the consistency of the numerical compu- 
tation. Then, let us redefine the fields in such a way that they will be approxi- 
mately of the same order of magnitude (this will avoid stiffness problems in the 
integration): 



where A, $ and £ will be the fields to integrate. It is also useful to set k = 
K/mpi and e c = e' c m pl . Once the equations have been rewritten under these 
considerations, we choose the initial conditions as follows: is set to 1 and 
4>i to zero (we therefore start with a "6are" gauge coupling constant equal to 
unity in the Einstein frame). We choose the ratio pBi{a — a,i)/e c = r so that we 
can control the type of gauge dynamics (BI, YM or transition) we start with. 
Then, we choose the value of (j)'(ai) — <j)'i so that the initial expansion rate will 
be given by Hf = n/3pBi(ai)/ (1 — n^^/Q). This gives also fa as it is equal to 
Hify. Then, without loss of generality, we can assume Ci = and determine &i 
from the postulated value of ps/(a — a{). Numerical integration of the system 
of six ODE's is performed using the standard method of Shampine-Gordon |27| . 
To monitor the accuracy of the numerical solution, we compute the absolute 
violation of the hamiltonian constraint (|10|l : 
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all along the integration. The numerical integration makes the violation of the 
hamiltonian constraint diverging exponentially with time and we indicated in 
the previous figures the final absolute error reached for each of the numerical 
solution that were presented. 
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